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In the paper, we have proposed and explored a way of creating a Rashba-Dresselhaus spin- 
orbit coupling (SOC) in an atomic-molecular Bose-Einstein condensates (BECs) by combining 
the spin dependent photoassociation and Raman coupling. Separately, the Raman coupling, which 
induces SOC 1 and photoassociation, which produces molecular BEC 7 ®, are well established 
experimental techniques, but there was no previous theoretical or experimental attempt to combine 
them. Amongst other results, we find a new type of topological excitation, the carbon-dioxide-like 
Skyrmion, in which two half-Skyrmions of the molecular BECs couple with one Skyrmion of the 
atomic BECs to form a carbon-dioxide-like structure. Our results show that the SOC can control 
the formation and distribution of the carbon-dioxide-like Skyrmion. This kind of Skyrmion can be 
measured via the vortices structures using the time-of-flight absorption imaging technique in the 
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real experiment. 


In this Supplementary Information, we first show the effects of atom-molecule conversion 
on the SOC atomic-molecular BECs, including consideration of the densities and the positions of 
vortices. Furthermore, we explore the effect of Raman coupling on the spin textures in the SOC 
atomic-molecular BECs. We prove that the carbon-dioxide-like Skyrmion would not occur in the 
normal rotating spinor atomic-molecular BECs without referring to both the Raman coupling and 
SOC. Additionally, we demonstrate that the transformation (S^., S^, S^) = (S~, S y , —S x ) does not 
change the topological charge. Finally, the solutions for characterizing the carbon-dioxide-like 
Skyrmion are provided. 

The effect of atom-molecule conversion on the vortices in the spin-orbit coupling atomic- 
molecular BECs. Now, we perform some experiments where the atom-molecule conversion is 
changed to be 0, 0.01 and 0.02. Figures Sla-Slc show the densities of the SOC atomic-molecular 
BECs which is obtained under the equilibrium state. The five components are pointed out by the 
title: | t), | {), m F = — 1 , m F = 0 and m F = +1, respectively. Similar to the experiment in 
Fig. 3, there are lots of vortices inducing. Figures Sld-Sll indicate the position of vortices. Here, 
the first column only shows the position of vortices formed by the atomic BECs (i.e., the spin-up 
and down components). The second column shows those formed by the molecular BECs (i.e., 
m F = 0, ±1 components). The third column indicates all the positions of vortices. Note that we 
use three rectangles (the red, blue and black) to demonstrate the relationship between vortices and 
that among molecular vortices. In all of these plots, the three different rectangles are of the same 
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sizes and the same positions, respectively. 


It is easy to find that each vortex of spin-up component is just above a coupled vortex of 
spin-down component in Fig. Sid. Off the x-axis, the distances between the two coupled vortices 
is very close (see the black rectangle), they become bigger in the blue rectangle, and they are the 
biggest in the red rectangle. These phenomena also occur in Figs. Slg and Sli. In Figs. Sle, 
Slh and Slk, the molecular vortices tend to form the three vortices structure where vortices of the 
m F = +1, rn F = 0, and m F = — 1 components couples and arranges in sequence. The black 
and blue rectangles in Figs. Sle and Slh do not catch some three vortices structures. But they can 
catch some three vortices structures in Fig. Slk. Comparing the results in Figs. Slf, Sli and S11, 
we can see clearly some vortex clusters and the vortices lattice (see the red rectangles) only in Fig. 
S11. Therefore, the vortices position indicates that the creation of vortices cluster depends on the 
atom-molecule conversion. 

The atom-molecule conversion enhances the coupling between atomic vortices and molecu¬ 
lar ones. When this conversion is terminated (i.e., x — 0-0), there is no coupling among atomic 
vortices and molecular vortices. In this situation, we obtain the regular distribution of atomic vor¬ 
tices in the black and blue rectangles and the relatively disordered molecular vortices in Fig. Sle. 
As the increase of the atom-molecule conversion, the molecular vortices have to unite with the 
atomic vortices and they are dragged into the black and blue rectangles. Thus, we can see the 
vortex clusters occurring in Fig. S11. This effect also happens in the red rectangles and leads to the 
creation of vortices lattice. Therefore, we can conclude that the atom-molecule conversion plays a 
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role in coupling the atomic vortices and the molecular vortices to form the vortex clusters and the 
vortices lattices. 

In real experiment, the photoassociation process can cause the atom-molecule conversion 
7 l3 . It is implies that the the photoassociation process are terminated if we choose the coeffi¬ 
cient of atom-molecule conversion to be 0. In that case, it is impossible to obtain any molecular 
BECs as well as the corresponding topological excitations in them. In this study, our method of 
damped projected Gross-Pitaevskii equations^ will behavior just likes to catch the particles with 
a shallow trap and exchange particles with an external thermal reservoir. This process is kept on 
until the system reaches the equilibrium state. So the molecular BECs is not completely coming 
from the trapped atomic BECs via the atom-molecule conversion. However, our method can well 
demonstrate the effect of atom-molecule conversion on the vortices and Skyrmions in the SOC 
atomic-molecular BECs, because the external thermal reservoir only provides the particles for the 
system but not brings in any other effects. Thus, the results of x — 0 can be viewed as an ex¬ 
periment where one firstly have gotten a mixed atomic-molecular BECs, then the lasers for the 
spin-dependent photoassociation are turned off and the atomic-molecular BECs are further cooled 
to the ground state. Therefore, the results are credible in Figs. Sla, Sid, Sle, Slf and in Figs. 5a, 
5d, 5g, 5j, 5m. 

The effect of Raman coupling on the carbon-dioxide-like Skyrmion in the spin-orbit coupling 
atomic-molecular BECs. To clarify why the carbon-dioxide-like Skyrmion can be induced in the 
atomic-molecular BECs, we further investigate the effect of Raman coupling Qa on the vortices 
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pattern and spin textures. In Figs. S2a-S2c, we select the strength of the Raman coupling with the 
value of 0, 0.2a; and 0.6a;, respectively. Comparing these results, the phase separation between the 
two components of the atomic BECs is the most obvious when the strength of the Raman coupling 
decreases to be 0. The reason is that Raman coupling can make the spin-up and spin-down atoms 
convert with each other. When we terminate this term = 0), the phase separation between 
two components of atomic BECs becomes the most obvious. Thus, there is little overlapping 
of densities between the two components of atomic BECs. Similarly, the phase separation also 
happens in the molecular BECs, where the mp = 0 component is in the middle. The corresponding 
vortices positions are shown in Figs. S2d-S2f. When the strength of the Raman coupling is 0 (see 
Fig. S2d), the carbon-dioxide vortices [the black rectangle and blue rectangle show two examples] 
are induced, and the regular structure of vortices lattices near the .x-axis disappears [i.e., there is no 
vortices in the red rectangle]. When the Raman coupling increases to be 0.2c o (see Fig. S2e), the 
red rectangle indicates that the regular structure of vortices lattices occurs near the x-axis. At the 
same time, vortices begin to form some clusters as shown by the black and blue rectangles. Note 
that not all the vortices form the clusters. When the Raman coupling is up to be 0.6a; (see Fig. S2f), 
the vortices form lattices in the red rectangle and the formation of vortex clusters becomes more 
clearly than that in Fig. S2e. Therefore, Raman coupling induces the coupling between atomic 
vortices of spin-up and spin-down components. Certainly, it also induces the coupling among the 
molecular vortices of the = 0, ±1 components. 

Raman coupling plays an important role in inducing the carbon-dioxide-like Skyrmion in 
atomic-molecular BECs. Figures S3a-S3c (Figures S3d-S3f) show the corresponding spin textures 
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of atomic BECs (molecular BECs) in Figs. S2a-S2c, respectively. The spin textures are obtained 
under the transformation (S^., S^. S^) = (S,.. S ;l/ . —S x ). We use the blue ellipses to mark out the 
carbon-dioxide-like Skyrmions. In Figs. S3a and S3d, there are circular structures of spin textures 
between the red lines. However, we would not like to view them as the effective Skyrmion struc¬ 
tures. Figures S4a and S4d show the topological charge densities of the spin textures in atomic 
BECs and those in molecular BECs. The topological charge density in Fig. S4a shows that the 
Skyrmion structure in atomic BECs off x-axis is not very clear, though the carbon dioxide vortices 
exist. The main reason is that it is difficult for atomic BEC of spin up (down) components filling 
in the atomic vortices of spin down (up) component in the absence of Raman coupling. While 
the spin textures along x-axis come from the domain wall of the atomic BECs. The topological 
charge density in Fig. S4a demonstrates a channel with a negligible change along the x-axis. In 
fact, there is little vortices structure but the domain wall. Here, we would not like to discuss this 
special structure. Analogous cases also occur in Figs. S4d and S4g. Certainly, the half-Skyrmions 
occur in Fig. S4d off x-axis. Figure S4g further proves that the structure of the carbon-dioxide¬ 
like Skyrmions is not very clear. All these results are caused by the absence of Raman coupling. 
It is difficult for the components of atomic BECs (molecular BECs) to fill with each other. Thus, 
normal Skyrmions are restrained though some atomic vortices are induced. Naturally, it is hard to 
obtain the clear carbon-dioxide-like Skyrmions. 

Raman coupling improves the creation of carbon-dioxide-like Skyrmions. In Fig. S3b, the 
structure of the Skyrmion (see the blue ellipses) is much more clear than that in Fig. S3a. We 
also use the red lines to indicate the region along the x-axis. The topological charge densities in 
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Fig. S4b show that there are Skyrmions occurring along between the two red lines. However, at 
the same region we can not find the singular points in Fig. S4e. Thus, there is no complete half- 
Skyrmion induced along the x-axis. When the strength of Raman coupling increases to 0.6a; (see 
Figs. S3c, S3f, S4c, S4f and S4i), there are complete half-Skyrmions whose topological charge 
densities have a singular point along the x-axis. Therefore, we can use the blue ellipses to mark 
out all the structures as the effective carbon-dioxide-like Skyrmions. These results prove that the 
Raman coupling enhances the formation of carbon-dioxide-like Skyrmions. 

The ground state of the rotating atomic-molecular BECs without spin-orbit coupling. We 

now discuss the ground state of the rotating atomic-molecular BECs without SOC. In our exper¬ 
imental setup [see Fig. 1 in the paper], the SOC originates from the counter propagating Raman 
lasers. If we turn off the Raman lasers, the experimental setup would only induce the normal spinor 
atomic-molecular BECs. Our purpose of this section is to prove that it is impossible to induce the 
carbon-dioxide-like Skyrmion in the rotating atomic-molecular BECs without SOC, i.e., without 
the counter propagating Raman lasers. Differed from the case in Fig. 3a and Fig. SI a, both SOC 
and the Raman coupling are terminated. Figures S5a and S5b show the densities distribution and 
the corresponding phases, respectively. The vortices positions are indicated in Figs. S5c-S5e. The 
spin textures show that there is no half-Skyrmions in Figs. S5f and S5g. All the excitation of the 
spin textures are normal Skyrmions. In addition, the spin textures in atomic BECs are very similar 
to that in molecular BECs. The topological charge densities in Figs. S5h and S5i further prove that 
there is no half-Skyrmions. Meanwhile, the sum of the topological charge density does not show 
that there are any coupling between the spin textures of atomic BECs and molecular ones. There- 
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fore, the carbon-dioxide-like Skyrmion does not occur if the Raman lasers which induce SOC are 
not considered. 


Topological charge Q is unchanged under a transformation (S x , S y . SJ = (S z , S y , — S x ). We 
now demonstrate that the transformation (S x , S^, S 2 ) = (S z , S y , —S x ) does not change the value of 
the topological charge Q. It is well known that the topological charge density q is defined as 
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ds <9s, 


q = 4 


( 1 ) 


where s = S/|S|, and Q — f f qdxdy. For clarity, we use q(s x , s y , s,) to describe the normal 
expression of the topological charge q of Eq. ©, and q(s z ,s y , —s x ) to denote the topological 
charge under the transformation (S^, S' S^) = (S~, S^, —S x ). 


We can expand the expression q(s x , s y , 
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Thus, the transformation (S^., S y , S*J = (S z , S y , —S x ) does not change the value of the topological 
charge density q. Certainly, the topological charge 0 also is not changed under the transformation. 


Solution of carbon-dioxide-like Skyrmion. Here, we illuminate the solution of the carbon-dioxide¬ 
like Skyrmion. For clarify, we take the results in Fig. 2 as an example to discuss the solutions 
detailedly. We find that the Skyrmion of atomic BECs in Fig. 2 can be represented approximately 
by 


(0.75a;) 2 + y 2 - 4 x 0.2 2 e“ 22x [(o ^sa:) 2 ^^ 2 ]/ 2 
Lx “ “ (0.75a;) 2 + y 2 + 4 x 0.2 2 e - 22 >< [(o. 75 ^) 2 + y 2 ]/2 , 

4 x 0.2 x o.75xe- 22x [(°- 753; ) 2+ ^]/ 2 
ly ~ (0.75a;) 2 + y 2 + 4 x 0.2 2 e ~ 2 ^m^) 2 W ]/2 ’ 
-4 x 0.2 x y e -#AQ-T5x)*W\l* 

Slz = (0.75a;) 2 + y 2 + 4 x 0.2 2 e- 2ax K°- 75a 0 2 +J' a ]/ 2 ' 


Figure S6a plots the structure of the spin texture with Eqs. d2j). Although the structure is different 
from the normal circular Skyrmion, it can transfer to the circular one under the transformation 
(S^., S ', S^) = (S z , Sy, —S x ) [see Fig. S6b]. The topological charge density in Fig. S6c show that 
this spin texture are nontrivial Skyrmion with Q = — 1. In addition, we find that all the properties 
of this Skyrmion is very similar to that in Fig. 2. Thus, Eqs. © can well describe the center 
Skyrmion in the carbon-dioxide-like Skyrmion. 
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A single half-Skyrmion of molecular BECs in Fig. 2 can be represented by 


s 2y = _- exp[—2.5 2 (x 2 + y 2 )/2], 

v ■ x + y 

S 2 z = - / y exp[—2.5 2 (x 2 +y 2 )/2], 

V x 2 + y 2 

*2x = -^/i-sly-sL- 


( 3 ) 


Clearly, the two vectors s 2y , s 22 can vary from —1 to 1, but the vector s 2x varies from —1 to 0. 
This means that the spin vector can cover half of a unit sphere. By calculating the topological 
charge Q — f f s ■ (J| x j^)dxdy, where s comes from Eqs. ©, we find the topological charge 
approaches Q = —0.5. 


The two half-Skyrmion of the carbon-dioxide-like Skyrmion in Fig. 2 can be approximately 


written as 


S 3y — 
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V 
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(4) 


The structure of spin texture with Eqs. © is shown in Fig. S6d. If we perform the transformation: 
(S^., S ' , S).) = (S 2 , S y , — S x ), we can obtain two clear circular structure [see Fig. S6e]. Comparing 
it with Fig. 2f, the formation of the two spin texture is very similar. This indicates the half- 
Skyrmion can be well described by Eqs. © in our simulation. Furthermore, the topological 
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charge density in Fig. S6f is very similar to that in Fig. 2j. Therefore, Eqs. © can well describe 
the two half-Skyrmions in the carbon-dioxide-like Skyrmion. 
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Supplementary Figure SI The effect of atom-molecule conversion on the vortices in the 
SOC atomic-molecular BECs. The rotation frequency is Vt rotation = 0.5a;, where u = 200 x 
27 tHz. We set the parameters A = 1, f 1 A = 0.8, < 7 ^ = = g a with the scattering length 

a A = 101.8a B , gu = 0.95 g A , g n = 2g A , g s = 0.03g A , g A M = 0.5 g A , £ = 0 and = 2 /j a = 
3.5/r m , 0 = 7/i aj0 = 42fta;. (a-c) show the densities of each component when the system reaches 
the equilibrium state. The components are pointed out by the title: | t), | |), mp = —1, m F = 0 
and nn F = +1, respectively, (a) x — 0-0, (b) x — 0-01, and (c) x — 0.02. (d-l) The position of 
vortices. The title ‘atom’, ‘molecule’ and ‘total’ indicate that the first, second and third column plot 
the position of vortices in atomic BECs, molecular ones and the both, respectively, (d) The position 
of vortices formed by atomic BECs in (a), (e) The position of vortices formed by molecular BECs 
in (a), (f) shows all vortices’ position in (a). The red, blue and black rectangles are used to mark 
out that atomic vortices tend to march with each other, and distribute parallel with the x axis. 
Similarly, (g-i) show the corresponding position of vortices in (b). (j-1) show the corresponding 
position of vortices in (c). The particles number of each component (N^, N±, N_ u N 0 , N +1 ) are 
(6.7 x 10 3 , 7.7 x 10 3 ,8.0 x 10 3 ,3.8 x 10 3 ,8.1 x 10 3 ) in (a), (6.9 x 10 3 , 7.8 x 10 3 ,8.0 x 10 3 ,3.8 x 
10 3 , 8.1 x 10 3 ) in (b), and (7.0 x 10 3 , 7.8 x 10 3 , 7.9 x 10 3 ,4.0 x 10 3 , 7.9 x 10 3 ) in (c). The unit 
of length, strength of x, and SOC are 0.76 gm, frw, u, and 0.96 x 10~ 3 m/s, respectively. 

Supplementary Figure S2 The effect of Raman coupling on the vortices in the SOC 
atomic-molecular BECs. The rotation frequency is O rotation = 0.5a;, where oj = 200 x 27 tHz. 
We set the parameters A = 1, g^^ = g^^ = g A with the scattering length a A = 101.8a s , g^ = 
0.95 g A , g n = 2 g A , g s = 0.03^, g A M = 0.5 g A , x = 0.08, e = 0 and g m = 2/i a = 3.5/x mj0 = 
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7/Vo = 42/V. All the results are obtained when the system reaches the equilibrium state. The 
five components are pointed out by the title: | t), | |), rn F = — 1 , m F = 0 and rn F = +1, 
respectively, (a-c) The effect of Raman coupling on the densities of the SOC atomic-molecular 
BECs. (a) fl A = 0, (b) fl A = 0.2, and (c) fl A = 0.6. (d-f) The corresponding position of vortices 
in (a-c), respectively. The blue ellipses point out the regions of the carbon-dioxide-like Skyrmion. 
In addition, we mark the position of vortices in order to illuminate the relationship between spin 
texture and position of vortices clearly. The meanings of the marks are the same as that in Figs. 
S2d-S2f. The particles number of each component (A^, jVj., iV_i, N 0 , N + i) are (8.0 x 10 3 , 8.6 x 
10 3 , 8.6 x 10 3 , 2.1 x 10 3 , 8.6 x 10 3 ), (7.8 x 10 3 , 8.8 x 10 3 , 8.6 x 10 3 , 2.4 x 10 3 , 8.3 x 10 3 ), and 
(8.0 x 10 3 , 8.6 x 10 3 , 8.1 x 10 3 , 3.4 x 10 3 , 8.1 x 10 3 ), respectively. The unit of length, strength of 
X, fU and SOC are 0.76/im, hu, u>, and 0.96 x 10~ 3 m/s, respectively. 

Supplementary Figure S3 The effect of Raman coupling on the spin textures in the SOC 
atomic-molecular BECs. The spin textures are obtained under the transformation (S^, S' y , Sj = 
(S-, Sj,, —S x ). The color of each arrow indicates the magnitude of S' z . (a, d) The corresponding 
spin textures of atomic BECs and molecular BECs in Fig. S2a. (b, e) The spin textures of atomic 
BECs and molecular BECs in Fig. S2b. (c, f) The spin textures of atomic BECs and molecular 
BECs in Fig. S2c. The blue ellipses point out the regions of the carbon-dioxide-like Skyrmions. 
The position of vortices are marked out with the same symbols as that in Figs. S2d-S2f. The unit 
of length and SOC are 0.76/tm and 0.96 x 10~ 3 m/s, respectively. 

Supplementary Figure S4 The effect of Raman coupling on the topological charge den- 
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sities in the SOC atomic-molecular BECs. (a-f) The corresponding topological charge densities 
of the spin texture in Figs. S3a-S3f, respectively, (g) The total topological charge densities of (a) 
and (d). (h) The total topological charge density of (b) and (e). (i) The total topological charge 
density of (c) and (f). The blue ellipses point out the regions of the carbon-dioxide-like Skyrmions. 
In the region of x < 0, the position of vortices are marked out. The meanings of the marks are 
the same as that in Figs. S2d-S2f. The unit of length and SOC are 0.76/xm and 0.96 x 10 _3 m/s, 
respectively. 

Supplementary Figure S5 The ground state of rotating atomic-molecular BECs without 
SOC. The rotation frequency is O rotation = 0.5a;, where ui = 200 x 27 tHz. We set the parameters 
A = 0, tl A = 0, 0 tit = g u = g A with the scattering length a A = 101. 8a B , gu = 0.95 g A , 
g A M = 0.5 g A , g n = 2 g A , g s = 0.03 g A , x = 0.08, e = 0 and = 2g a = 3.5g mfi = 7g a . 0 = 
28 huj. (a) Densities of each components under the equilibrium state. Note that the components 
are pointed out by the title: j |), j |), m F = —1, m F = 0 and m F = +1, respectively, (b) The 
corresponding phases of each components in (a), (c-e) The position of vortices, (c) The position 
of vortices formed by atomic BECs. (d) The position of vortices formed by molecular BECs. 
(e) shows all vortices’ position, (f) The spin textures of atomic BECs under the transformation 
(S^., S ', S^) = (S z , S y, — Sz). (g) The corresponding spin textures of molecular BECs. The color 
of each arrow indicates the magnitude of S' z . (h) The topological charge densities of the spin 
texture in Figs. S5f. (i) The topological charge densities of the spin texture in Figs. S5g. (j) The 
total topological charge densities of (h) and (i). The particles number (A^, N_i, N 0 , N +1 ) is 
(2.72 x TO 3 , 3.35 x TO 3 , 2.48 x 10 3 , 2.80 x 10 3 , 2.23 x 10 3 ). The unit of length, strength of y, 
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and SOC are 0.76 fim, ftw, oj, and 0.96 x 10 3 m/s, respectively. 


Supplementary Figure S6 Solution of carbon-dioxide-like Skyrmion. (a) Structure of a 
Skyrmion with Eqs. ©. The color of each arrow indicates the magnitude of S\ z . (b) Structure 
of a Skyrmion with Eqs. © under the transformation (S^S^S'j = (S^S^, —S x ). The color 
of each arrow indicates the magnitude of —S± x . (c) The topological charge density q with Eqs. 
©. (d) Structure of two half-Skyrmions with Eqs. ©. The color of each arrow indicates the 
magnitude of s^ z . (e) Structure of two half-Skyrmions with Eqs. © under the transformation 
(S^, S ', s!j = (S 2 , S y , —S x ). The color of each arrow indicates the magnitude of — s 3x . (f) The 
topological charge density q with Eqs. ©. 
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Atomic-molecular Bose-Einstein condensates (BECs) offer brand new opportunities to rev¬ 
olutionize quantum gases and probe the variation of fundamental constants with unprece¬ 
dented sensitivity. The recent realization of spin-orbit coupling (SOC) in BECs provides a 
new platform for exploring completely new phenomena unrealizable elsewhere. However, 
there is no study of SOC atomic-molecular BECs so far. Here, we find a novel way of creat¬ 
ing a Rashba-Dresselhaus SOC in atomic-molecular BECs by combining the spin dependent 
photoassociation and Raman coupling, which can control the formation and distribution of a 
new type of topological excitation - carbon-dioxide-like Skyrmion. This Skyrmion is formed 
by two half-Skyrmions of molecular BECs coupling with one Skyrmion of atomic BECs, 
where the two half-Skyrmions locates at both sides of one Skyrmion, which can be detected 
by measuring the vortices structures using the time-of-flight absorption imaging technique 
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in real experiments. 


The production of atomic-molecular Bose-Einstein condensates (BECs) via the Raman pho¬ 
toassociation 17 has revolutionized the field of ultracold quantum gases and generated much inter¬ 
est in controlled molecular dynamics, ultrahigh-resolution molecular spectroscopy and quantum 
information processing, where molecules offer microscopic degrees of freedom beyond atomic 
gases. Atomic-molecular BECs possess unique properties that may allow for the study of new 
physical phenomena and lead to discoveries, reaching far beyond the atomic BEC. Furthermore, 
atomic-molecular BECs exhibit the non-Arrhenius superchemistry, which describes the coherent 
stimulation of chemical reactions via macroscopic occupation of a quantum state by a bosonic 
chemical species 1 10 . Recently, spin-orbit coupling (SOC) has been realized for ultracold gases 
of atomic bosons 11 16 and fermions^ 17 1S . Especially, SOC is ubiquitous in physical systems. This 
raises a question how to create SOC in the mixed atomic-molecular BECs. At least, the hyperfine 
states of ultracold molecules® 3 21 demonstrate that molecular systems possess the essential char¬ 
acteristics for the creation of SOC. However, there is no study of the SOC atomic-molecular BECs 
so far. 


SOC describes the interaction of the particle’s spin with its orbital motion. It represents 
a fascinating and fast-developing area of research significantly overlapping with traditional con¬ 
densed matter physics, but importantly containing also completely new phenomena unrealizable 
elsewhere. Novel effects have been reported, such as formation of the Bose mixture exhibiting a 
phase-separated propagation 21, as well as appearance of the plane wave and stripe phases^ 22 24 and 
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other nontrivial structures ^ 25 36 in the SOC atomic BECs. A number of schemes have been pro¬ 
posed to create general gauge fields and SOC for ultracold atoms G 7 5H. Further exploration of SOC 
on atomic-molecular BECs is, however, still lacking. Differing from the normal SOC ultracold 
atoms, the SOC atomic-molecular BECs would initiate the exploration of SOC to the untouched 
hybrid quantum systems. 

Here, we find a novel way of creating a Rashba-Dresselhaus SOC in the atomic-molecular 
BECs by combining the spin dependent photoassociation 21 42 and Raman coupling 1 By ap¬ 
plying this SOC, we demonstrate that the atomic-molecular BECs support a new type of topo¬ 
logical spin texture - carbon-dioxide-like Skyrmion, which is formed by two half-Skyrmions of 
the molecular BECs coupling one Skyrmion of the atomic BECs, where the former two half- 
Skyrmions locates at both sides of the Skyrmion. Our results not only indicate the immanent 
mechanism of a single carbon-dioxide-like Skyrmion, but also show that SOC can control the 
carbon-dioxide-like Skyrmions to distribute in parallel along the x axis. The SOC even can de¬ 
stroy the half-Skyrmion in the carbon-dioxide-like Skyrmions which are in parallel along the x axis 
when its strength is very strong. In addition, our results prove that the atom-molecule conversion, 
which describes the superchemistry, causes the combination of the Skyrmion and half-Skyrmions 
in this carbon-dioxide-like Skyrmion. If there is no the SOC or the accompanying Raman cou¬ 
pling, the carbon-dioxide-like Skyrmion will not occur. Furthermore, we find that the formation of 
the carbon-dioxide-like Skyrmion relates to some vortices structures. For instance, the Skyrmion 
in the carbon-dioxide-like Skyrmion may originate from a dipole structure of a pair of vortices in 
the atomic BECs, and the half-Skyrmion in the carbon-dioxide-like Skyrmion may derive from a 
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three-vortex structure, where two vortices of m F = —1 and rn F — 1 components are respectively 
on both sides of a vortex of m p = 0 component in spin-1 molecular BECs. Therefore, the detection 
of carbon-dioxide-like Skyrmion can be based on the vortices structures by using the absorption 
imaging technique in a real experiment. 

Results 

Creation of Rashba-Dresselhaus spin-orbit coupling in atomic-molecular BECs. Experimen¬ 
tally, the Rashba-Dresselhaus SOC can be obtained by using the Raman-coupling of two hyperfine 
states in the ground state manifold of alkali atoms 2H03 In that case, the strength of the SOC 
depends on the wave vectors of the Raman lasers. On the other hand, the effective out-of-plane 
(in-plane) Zeeman field is proportional to the Rabi frequency of the Raman process, i.e. to the ma¬ 
trix element of the two-photon coupling between a selected pair of atomic magnetic sublevels of 
the ground-state manifold. It is well known that the molecular BEC can be created by the Raman 
photoassociation of atoms in a condensate 1 7 . With the spin dependent photoassociation 21 42; , the 
hyperfine states of molecules are related to those of the converted atoms. Furthermore, a number 
of experiments have demonstrated a high degree of control over the hyperfine state of ultracold 
molecules via two-photon optical Raman transitions 21 42 ®3. Here, when considering the SOC in 
atomic-molecular BECs, we are dealing with the atoms characterized by two different internal 
(quasi-spin) states, as well as molecules with three spin states (F = 1, m F = 0, ±1). 

We consider a Bose gas of ultracold s ' Rb atoms, commonly used in both experiments on the 
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SOC atomic BEC 11 ® and the atomic-molecular BECs^E, to demonstrate the scheme for creating 
SOC atomic-molecular BECs [see Fig. la]. We take a BEC of up to 3 x 10 5 atoms, and select 
two atomic internal spin states with rrif = —1 and rrif = 0 of the s 'Rb 5S1/2 ground state / = 1 
manifold, to be referred to as the (pseudo-) spin up (| t)) and down (| ],)) states [see Fig. lb]. Here, 
/ and nrif are the quantum numbers of the atomic total spin and its projection along a quantisation 
axis. A pair of counter-propagating Raman lasers ( ujl 3 and c Ola) couple these spin states with the 
matrix element (Raman Rabi frequency) inducing the SOC of the atomic BECs [see Fig. lb]. 
The third spin state with rrif = 1 is detuned away from the Raman resonance due to the quadratic 
Zeeman shift and hence can be neglected. On the other hand, the counter-propagating Raman 
lasers ( u>l 5 , c olg and c olt) are used to induce the SOC for the spin-1 (F = 1) molecular BECs [see 
Fig. lc]. An extra laser ul 7 is needed to couple the molecular states with mj? = 0 and m F = 1 
detuned from the resonance between the states = -1 and rri F = 0 due to a quadratic Zeeman 
effect. To prevent population of the excited state the laser fields are to be sufficiently detuned from 
the single photon resonance. The detailed processes for characterizing the SOC atomic-molecular 
BECs are indicated in the Methods. 

The molecular BECs of s 'Rb 2 with hyperfine state |F = 1, m F = 0, ±1) are produced by 
the spin dependent photoassociation of atoms ^ 21 42 (see Figs. Id-If). A pair of atoms with various 
combinations of internal states is coupled to an excited molecule state | (p em ) via a laser beam ujli- 
Here, the photoassociation through the total spin-2 scattering channel gives molecular F numbers 
of 2. The three spin dependent photoassociation process can be written as: 11, 0)<g)|l, 0) —)• |2,0), 

11, —1) £5 11, 0) —> 12, —1) and 11, —1) ® 11, —1) —» |2, —2). Under the transition A F = —1 
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and A m F = +1, the excited molecule states further change into \<f> em ) with F — 1 and / = 3 Eh. 
where F, mp and I denote the total spin, its projection and nuclear spin quantum numbers of a 
molecule. Finally, \<j) em ) is coupled to a ground molecule state \(j) m ) via another laser beam uj L2 
[see Figs. Id-If]. As a result, the two-component atomic BEC can be coherently converted to the 
three-component (spin-1) molecular BEC with three internal spin states |F = 1, nip = 0 , ±1). 


The SOC atomic-molecular BECs can be described by (see Methods): 
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xVtVt+^o, and J\_i = -^^| + (£+M J 4A 2 )0_i. Here a =|, i and m = 0, ±1 labels, respectively, 
the atomic and molecular spin, d = 1,2,3 refers to the dimension of the system, M A ( Mm ) is the 


mass of the atom (molecule), F a=XtVjZ ( a a = x ,y,z ) are the spin-1 (Pauli) matrices, g a ^ q = 
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are, respectively, the strengths of atom-atom and atom-molecule 


MaMm/(Ma+Mm) 

interaction, g n and g s are the strengths of the spin-independent and spin-dependent interactions 
between the molecules ®2EID. Here also ip a and 0 m denote the macroscopic wave functions of the 
atomic BEC in the internal state |a) (a =j\ |) and those of the molecular BEC in the internal state 
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| m F = m) (m = 0, ±1). The parameters A, Q-i, x, e, and c o are the strengths of SOC, Raman 
coupling, atom-molecule conversion, Raman detuning, and trap frequency, respectively. 

A carbon-dioxide-like Skyrmion in the SOC atomic-molecular BECs. Spin texture excitations, 
such as Skyrmion and half-Skyrmion (meron), are particle-like topological entities in continuous 
field and are recognized to play an important role in many condensed matter systems ^ 47 52 . The 
SOC atomic-molecular BECs provide a new platform for exploring the spin texture excitations 
coming from atomic BECs 55 and molecular BECs simultaneously. Now, we firstly present 
a way to obtain a single carbon-dioxide-like Skyrmion, which consists of one Skyrmion with a 
topological charge —1 in the atomic BECs and two half-Skyrmions with a topological charge — | 
in the molecular BECs. On both sides of the Skyrmion are the two half-Skyrmions and they couple 
as a unity resembling the carbon dioxide molecule. 

In numerical simulations, we set the initial chemical potential // rr ,.o = 2/x a ,o = 8ftw. The 
final chemical potential of the noncondensate band is increased to /i m = 2/i a = 28hu. We use 
the parameters of atomic-molecular BECs of 87 Rb with M m = 2M A (M,\ = 144.42 x 10 2 ' Kg), 
the trapping frequency ui = 200 x 27 tHz, and strength of atom-atom interaction g^ = g A with 
the scattering length a A = 101.8ob, where aj> is the Bohr radius. In addition, when the change 
in energy of converting two atoms into one molecule (AC/ = 2 Uto, ~ Urm) *1 without including 
internal energy, approaches zero, we can obtain the value 2 gA = 9 m, and we set g n = g M = 
2 g A . For simplify, we set the atom-molecule interaction g^-i = (?t,o — fl't.+i = 9i- 1 = 9i, o — 
gi t+ i = gAM■ The unit of length, time, and energy correspond to a/ h/(M A u>) (~ 0.76/im), 
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(« 0.96 x 10“ 3 s), and hw, respectively. In the rotating frame, H is replaced by H — Vt rotation L z , 
where L z = — ih(xd y — yd x ) is the z component of the orbital angular momentum and £1rotation is 
the angular frequency of rotation. 

The densities and the corresponding phases of the SOC atomic-molecular BECs are obtained 
for the equilibrium state (see Figs. 2a and 2b). One can see that both components of atomic BECs 
contain a vortex excitation. The spin-up (spin-down) component fills in the vortex of the spin- 
down (spin-up) component. Thus, this structure is similar to a vortices molecule 56 . At the same 
time, two vortices are induced in each molecular component, respectively. The densities of all the 
components are symmetric with respect to the y axis. In addition, only the density of the m F = 0 
component is symmetric with respect to the x axis. In this way, these vortices form a cluster. 

The emergence of the vortices cluster results from a combination of the SOC, Raman cou¬ 
pling, atom-molecule conversion and rotation effect. The Rashba-Dresselhaus SOC causes a trans¬ 
lational displacement along the y axis on the rotating BECs. The translational displacement for 
different spin state is opposite. Thus, the spin-up and m F = +1 components trend to move to the 
y < 0 region, and the spin-down and m F = — 1 components go to the y > 0 region. Note that the 
m F = 0 component is not affected by the Rashba-Dresselhaus SOC. Rotation also induces vor¬ 
tices. In the atomic BECs, the spin-up and the spin-down component have to fill with each other 
in the region of the vortices. Similar case occurs in the molecular BECs with the vortices of the 
m F = 0 component locating at the center. The atom-molecule conversion forces the atomic vortex 
of spin-up (spin-down) component to couple with two molecular ones of mp = +1 (m F = —1) 
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components, just like that in the normal atomic-molecular BECs^ 7 . This leads to the formation of 
the above vortices cluster. 

Since our system includes the pseudo-spin-1/2 atomic BECs and the spin-1 molecular BECs, 
there will be spin textures coming from atomic BECs and molecular BECs, respectively. Thus, 
this mixed system will show the interaction of the spin textures induced by different hyperfine spin 
states. In pseudo-spin-1/2 BECs and spin-1 BECs, the spin texture 56 58 ^0i s defined by 

V , m(^)m,nV , n/|'0| 2 , (« = X,y,z), 

m,n=i'-i- 

sr‘“" te = E (3) 

m,n=0,d=l 

where cr Q and F a are the Pauli and spin-1 matrices, respectively. 

The spin textures in the atomic BECs and molecular BECs are shown in Figs. 2c and 2d, 
respectively. In Fig. 2c, the arrows form a circle located near the position (0, 0). The arrows trend 
to point down in the region y < 0 and up in the region y > 0. Similarly, in Fig. 2d, the arrows 
form two half circles located near the positions (0, —0.7) and (0, 0.7), respectively. Performing 
a transformation (S^., S^, = (S z , S y , —S x ), one can identify clearly a Skyrmion [see Fig. 2e] 

and two half-Skyrmions [see Fig. 2f]. In the Supplementary Information, we have proved that 
the transformation does not affect the value of the topological charge Q [Q = -^ J J q(x, y)dxdy, 
where q(x,y ) = s ■ (J| x ||), s = S/|S|]. Additionally, we provide approximate solutions for 
the carbon-dioxide-like Skyrmion [see the Supplementary Information]. Figure 2g is a combined 
scheme of a carbon-dioxide-like Skyrmion, where a Skyrmion of atomic BECs couples with two 
half-Skyrmions of molecular BECs. Figure 2h shows the differences of spin vectors between 
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the atomic BECs and the molecular BECs (y = 0). It indicates that the spin vectors of atomic 
BECs and molecular BECs would trend to be parallel to each other except at the region of the 
carbon-dioxide-like Skyrmion. Figures 2i and 2j show the topological charge density q(x. y). The 
positions of vortices are pointed out to highlight the relationship between Skyrmion and vortices. 
We find that the topological charge approaches —1 and — A for the two types of spin texture, 
respectively. The two zero points in the Fig. 2j further prove that we obtain two half-Skyrmions. 
Figure 2k shows the sum of topological charge density [ q sum (x , y) = q atom {x , y) + q molecule ( x , y)], 
which further clearly indicates the carbon-dioxide-like structure. 

The effect of spin-orbit coupling on carbon-dioxide-like Skyrmions in atomic-molecular BECs. 

We now explore an effect of the SOC on the rotating atomic-molecular BECs. Figures 3a, 3c, 3e 
display densities of the atomic-molecular BECs obtained for the equilibrium state at VL rota tion = 
0.5cu. Figures 3b, 3d, 3f demonstrate the corresponding phases. Without the SOC (A = 0), two 
atomic and three molecular components of the BECs do not display the phase separation [see Fig. 

3 a]. Yet there are a lot of vortices in each component, and they can be identified by the phase 
image [see Fig. 3b]. However, adding the SOC (A = 1), the two components of atomic BECs split 
into two parts [see Fig. 3c]. The component described by is located below the x axis, whereas 
-04, is localised mainly above it. Similarly, the phase separation also occurs in the molecular BECs. 
Because the spin-up (spin-down) atoms can form the 0 +1 (0_i) component, the corresponding 
molecular BECs of 0 + i (0_i) also tends to locate below (upper) y = 0 axis. The 0 O component is 
spread around the x direction and squeezed in the y direction, and the area is filled with the vor¬ 
tices. When the strength of the SOC is up to A = 2, the phase separation becomes more obvious 
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[see Fig. 3e]. In fact, the phase separation does not take place at all if the rotation is absent. There¬ 
fore, the combination of SOC and rotation induces the phase separation in the atomic-molecular 
BECs. 

To detect the underlying relationship between the topological excitations, Figures 3g-3i show 
the position of vortices, where the results correspond to Figs. 3a, 3c, and 3e, respectively. Note 
that we do not indicate the vortices where the densities of BECs are very low. In that case, the 
contribution of energy and momentum to the system almost can be neglected. Without the SOC 
[see Fig. 3g] five types of vortices form certain clusters (a black rectangle shows an example). 
The repulsion among vortices of each component forces them to spread over the BECs, but they 
do not form a regular lattice. At the same time, the repulsion among vortex clusters also leads to 
the deformation of the clusters shown in Fig. 3g. In Fig. 3h (A = 1), a black rectangle shows 
some regular structure of vortices occurring around the region along the x-axis. In fact, they are 
essentially vortex clusters which degenerate into a line. The number of vortices in each column is 
the same as that in Fig. 2a. The SOC leads the vortices to occur along the x axis, but the repulsion 
among vortices squeezes them to distribute in a line. Thus, several lines of vortices form a 2D 
lattice. In a periphery of the lattice there are the remaining vortex clusters. Outside the BECs, the 
vortices fail to form the full vortex clusters, because the phase separation causes some components 
to disappear locally. When the strength of the SOC is increased up to A = 2, the structure of the 
vortices lattice (a black rectangle region in Fig. 3i) is different from that in Fig. 3h. Here, each 
vortex of 0 O component tends to overlap with an atomic vortex. And the vortex number of (j > + 1 

component decreases to 1 in each column. The blue rectangle and the red rectangle in Fig. 
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3i indicate that vortices form some carbon dioxide structures further away form the center of the 
vortex lattice (black rectangle). Above the black rectangle, vortices in components of and ©_ | 
form the carbon dioxide structure, whereas the atomic vortex are situated in the middle part, and 
two molecular vortices are located at the two sides symmetrically. Similarly, the carbon dioxide 
structure of vortices in components of ^ and e> + \ occurs below the black rectangle. These results 
indicate that the strong SOC deepens the phase separation. At the same time, it forces the vortices 
cluster to transform into the lattice around the x axis and even into the carbon dioxide structures in 
the remaining regions. 

We now discuss the effect of the SOC on the spin texture on the above results. Figure 4 
shows the spin texture after the transformation (S^., S y , S J = (S 2 , S y , —S x ). In the Supplementary 
Information, it has been proved that the transformation only changes the formation of the spin 
texture but does not change its character. Figures 4a and 4d show the spin texture in atomic BECs 
and molecular BECs without SOC (A = 0), respectively. The topological charge density in Fig. 
4g (Fig. 4j) shows that they are Skyrmion (half-Skyrmion) with Q = — 1 (Q — — |). The position 
of vortices proves that the Skyrmion (half-Skyrmion) originates from the vortices dipole of atomic 
BECs (the three vortices structure of molecular BECs). The sum of topological charge density in 
Fig. 4m indicates that the Skyrmion of atomic BECs and the half-Skyrmion of molecular BECs 
form the carbon dioxide structure, which is marked out by blue ellipses. 

The SOC can redistribute the carbon-dioxide-like Skyrmions in the atomic-molecular BECs. 
For A = 1, we also obtain some Skyrmions and half-Skyrmions in the atomic BECs and molecular 
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BECs, respectively. Due to the phase separation caused by the SOC and rotation, some Skyrmions 
are forced to distribute along the x-axis in Fig. 4b. Similarly, some half-Skyrmions distribute 
parallel along the x-axis in Fig. 4e. The topological charge density in Figs. 4h, 4k and 4n shows 
the carbon-dioxide-like Skyrmions distribute in a line along the x-axis. 

The strong SOC can even redesign the spin excitations of the textures. When the strength 
of the SOC A reaches 2, the Skyrmion can occur in the atomic BECs [see Figs. 4c and 4i]. The 
Skyrmion along the x-axis is induced by the vortices dipole but away from the x-axis they are 
induced by a single vortex. In the molecular BECs, the half-Skyrmion is induced by a single 
vortex of </>_i component in the region of y > 1.5 [see Figs. 4f and 41]. Along the x axis (between 
the black lines), the spin texture transforms into a structure with a fractional topological charge. 
This structure is not constructed by the normal half-Skyrmion because the singular points in the 
topological charge densities disappear. The sum of the topological charge density in Fig. 4o 
further proves that the carbon-dioxide-like Skyrmion indeed occur in the region where y > 1.5 or 
y < -1.5. 

The effect of atom-molecule conversion on the carbon-dioxide-like Skyrmion in the spin-orbit 
coupling atomic-molecular BECs. The effect of atom-molecule conversion on the spin textures 
in atomic-molecular BECs is shown in Figs. 5a-5f. We find that it is hard to mark out the carbon- 
dioxide-like Skyrmion after the atom-molecule conversion being terminated (i.e., x — 0.0, see 
Figs. 5a and 5d). The topological charge densities in Figs. 5g and 5j prove that lots of Skyrmions 
and half-Skyrmions are induced. However, without the atom-molecule conversion, the two ki nds 
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of spin textures do not couple with each other. The sum of the two topological charge densities in 
Fig. 5m further shows that it is hard to separate these structures into independent unities, which 
can be marked out by the blue ellipses. When the atom-molecule conversion is x — 0.01, we 
can divide generally the topological charge densities into several independent unities by the blue 
ellipses (see Figs. 5h, 5k and 5n). This kind of dividing also is shown in Figs. 5b and 5e. Clearly, 
these structures are carbon-dioxide-like Skyrmions. In fact, we find some ellipses would overlap 
partly. When the atom-molecule conversion is up to x = 0.02, there is no any overlapping among 
the dividing (see Figs. 5c, 5f, 5i, 51 and 5o). Thus, the atom-molecule conversion, which is the 
key point of the superchemistry effect, induces the coupling between Skyrmion in atomic BECs 
and half-Skyrmion in molecular BECs. Therefore, the carbon-dioxide-like Skyrmion is a special 
phenomenon in atomic-molecular BECs because the atom-molecule conversion is a needed factor 
to induce the coupling. 

In the Supplementary Information, we indicate the corresponding densities distribution of 
the atomic-molecular BECs [see Figs. Sla-Slc], and the positions of vortices [see Figs. Sld-Sll]. 

Discussion 

We have explored the formation of the carbon-dioxide-like Skyrmions, where one Skyrmion of 
atomic BECs couples with two half-Skyrmions coming from molecular BECs to form the carbon- 
dioxide structure in the SOC atomic-molecular BECs. Our results also show systematically the 
effect of SOC and atom-molecule conversion on the creation of this spin textures excitation. This 


14 


kind of Skyrmions originate from some vortices structures such as the vortex cluster and the 
carbon-dioxide atomic-molecular vortices. SOC can make the vortex clusters favor to distribute 
along the x axis, and the repulsion among vortex clusters squeezes them into lines. Thus, the 
vortices lattice tends to occur around the x axis. Meanwhile, the carbon-dioxide-like Skyrmions 
are forced to distribute in parallel along the x axis due to the effect of SOC on the vortices. Very 
strong SOC changes the structure of vortices lattice and destroys the half-Skyrmion around the 
x axis. But it also makes the vortices cluster degenerate into carbon dioxide vortices, which can 
induce the carbon-dioxide-like Skyrmion. Skyrmions are often viewed as the ideal information 
carriers of great prospect in future application via the creation and annihilation 49 —. The carbon- 
dioxide-like Skyrmion is the result of combining the separated topological excitations to form a 
structural cell. This property implies that if the carbon-dioxide-like Skyrmions can serve as in¬ 
formation carriers 50 , including the creation and annihilation, their distribution in lines and in 
parallels, and the combination and separation will demonstrate the potential for information stor¬ 
age. Thus, the carbon-dioxide-like Skyrmion allows more complex control than normal Skyrmion 
discussed previously for information carriers. 

In the Supplementary Information, we further prove that the Raman coupling enhances the 
creation of carbon-dioxide-like Skyrmion (see Supplementary Figs. S2, S3 and S4). The case 
in Fig. 3(a) of single Raman coupling without SOC would be realized by lasers with the same 
direction. It is impossible to only hold the Raman coupling but neglect the SOC with the counter- 
propagating Raman lasers. If the counter-propagating Raman lasers are turned off, i.e., both Raman 
coupling and SOC are terminated, the carbon-dioxide-like Skyrmions do not occur in the rotating 
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atomic-molecular BECs although the superchemistry effect is considered (see Supplementary Figs. 
S5). Furthermore, there are no nontrivial spin texture excitations if the rotation frequency decreases 
to be zero. Due to the dependence of the spin-dependent photoassociation and Raman lasers, 
the carbon-dioxide-like Skyrmion is a special phenomenon existing in the SOC atomic-molecular 
BECs with hyperfine spin states. 

In real experiment, the carbon-dioxide-like Skyrmion can be identified by a direct observa¬ 
tion of the densities via time-of-flight absorption imaging technique®® 3 ; which can measure the 
characteristics of the carbon-dioxide-like Skyrmion by means of corresponding quantum vortices. 
The experiments should begin with nearly pure 87 Rb BEC of approximately 3.0 x 10 5 or more 
atoms in the \f = 1, my = —1) state 1 '-®. An initially off-resonant radio-frequency magnetic field 
can be used to prepare equal mixtures of \f = l,m/ = —1) and |/ = l,m/ = 0). Then, the 
spin-dependent photoassociation of atom 21 42 creates the molecules of |F = 2, rnp = 0, —1, —2), 
which are used to produce the excited molecules of |F = 1 ,mp = 0, ±1) via the transition 
A F = —1 and Amp = +1, and finally the excited molecules degenerate to be the molecular 
BECs of |F = 1 ,mp = 0, ±1) [see Figs. Id-If]. At the same time, the counter-propagating 
Raman lasers L3 and LA (L5, L6 and L7) shown in Fig. 1 are used to induce SOC for atomic 
(molecular) BECs^ 11 ®. Combining rotation effect and the further cooling process, one would ob¬ 
tain the SOC atomic-molecular BECs at the equilibrium state. In this study, our results have shown 
that the carbon-dioxide-like Skyrmion is related either to vortices cluster or to carbon-dioxide vor¬ 
tices. The vortex cluster consists of two atomic vortices and six molecular vortices [see the black 
rectangle in Fig. 3g]. The structure of the carbon-dioxide vortices is simple [see the red and blue 
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rectangles in Fig. 3i]. Therefore, one can apply a ballistic expansion (approximate 20ms) of the 
BECs in a gradient magnetic field to separate the different spin states Vortices will be im¬ 
printed by ramping the magnetic field B z —> 0. After the structures of vortices are measured, the 
carbon-dioxide-like Skyrmion can be identified. 

In summary, we have constructed an experimental proposal for creating SOC in atomic- 
molecular BECs by combining the spin dependent photoassociation and Raman coupling. This 
type of SOC will create the novel topological excitation - carbon-dioxide-like Skyrmion and also 
control the structures of this Skyrmions, which can be measured by the time-of-flight absorption 
imaging technique. This work opens a new window to reali z e of SOC in multi-component ultracold 
gases and manipulate complex nontrivial topological excitations in future experiments. 

Methods 

The equations of the spin-orbit coupling atomic-molecular BECs. The SOC of atomic BECs 
can be generated by two additional Raman beams ( uj L3 and which couple the atomic internal 
states | t) an d | |). The single-atom Hamiltonian induced by the Raman coupling is t)(l 
| e - 2 ik x x } l c ^ w j t j 1 () ( th e atomic Raman coupling intensity. Furthermore, Raman beams (u>l 5 , 
oj fjj and uj L7 ) couple different molecular internal states rn F = 0, ±1), and lead to Hamiltonain 
(| m F = 0 )(m F = —1| + | m F = 1 )(m F = 0|) fl M e lkxX + h.c., with X the position of the mass 
center of the two atoms and Vt M = v^Xt. 

To clearly illustrate effects of the SOC, we introduce a unitary transformation Ua = e lr7zkxX 


17 


and Um = ^ m=0 ±1 e 2zmkxX \rriF = rn)(rriF = m\ for the atomic states and molecular states, 
respectively, with a z = | t) (t | — I i) (I I • Therefore, in the spin rotated frame induced by a 
combined transformation U = UaP a + UmP m , the Hamiltonian is H = UII'U\ where IP is the 
Hamiltonian in the original (un-rotated) frame, P A and P M being the projection operators onto the 
atomic and molecular states. The straightforward calculation shows that, for the condensate in the 
d-dimensional space (d = 1,2, 3), in the second quantization we have 

H = Haf + HmF + Hsoc + PI Raman + PAM + Pi- (4) 

Here Haf = E a =t ,i I ( 2 M a) + M A u 2 r 2 /2}^ a (r)dr and H MF = £ m =o,±i J 4>l 

(f)[— K 2 W 2 d /(AM a) + Mau 2 t 2 + £ + MA\ 2 Tn 2 ](f) m (r)df are the free Hamiltonians for atoms and 
molecules, respectively, where Ma is the single-atom mass, u> is the trapping frequency, e is the 
difference between the binding energy to the two-atom molecule and the frequency difference of 
the two laser beams which induce atom-molecule coupling, the factor A is defined as A = Ic x /Ma, 
the operator ^ a (r) (a =|, |) is the creation operator for a single atom at position r with internal 
state |a), and ©J n (F) (m = 0, ±1) is the creation operator for a single molecular at position r with 
internal state | mp — m). In Eq. © the SOC is described by 

H S oc — A J ty{r){-itu) x )'ih{r)dr- A j {r){-ihd x )i)^{r)df 
+A J 4>\(r)(—ihd x )4>i(r)dr — X J ihd x )4>_i(r)dr. 

The Ramman coupling between internal states is described by 

H R aman = ^A J ^\{r)^{r)df ■ + Vt M J <j>\ (f)}o(r)df+ 0 M J (r)df+h.c.. 
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The laser-induced atom-molecule transition is described by 

H am = X fd^mV/rt + (r)0f(^) + 0-i /^ + h -c.]. 

In our system the inter-particle interaction is described by Hj in Eq. ©, and we have IIj = 
Uaa+Uam+Umm ■ Here U A a = E a =t,4 .9aa f dr^l{r) 2 ^ a (r) 2 /2+g n J dr^|(r)^j(r)^(f^^(r) 
describes the atom-atom collision, U AM = E Q =1U Em=o,±i 9 am f de¬ 

scribes the atom-molecule collision, and Umm describes the molecule-molecule collision. In usual 
case the collision between molecules with different internal states are quite complicated. Here 
for simplicity we assume in our system the scattering of three-component molecules has the same 
property as the scattering of F — 1 atoms & 5 E3, i.e., the scattering amplitude is totally determined 
by the total spin of the two molecules. In that case we have U M m = E m m'= o ±1 / (e) 

<i>m'{r)4>m{r)\ + Ec £m,n,fc,Z= Oil f {F a )nm (F a ) 0/ (r )}.. 


In this way, we obtain the coupled equations of the atomic-molecular BECs with SOC. It can 
be written as 

, h 2 X7j , M A LU 2 r 2 n 


in a Jn 

/it 


2 M a 


+ 


-]0f + t + 


0f,a|^a| 2 + ^t,m|0m| 2 )^t+ (^XV’t^’+l + XV’l^o), (5) 

a=t,4. m=0,±l 




M A w 2 r 2 


2M a 




+ (J]^“l^ a | 2+ S ^4.,m|^m| 2 )^4.+ (V^X^-l + X^o), ( 6 ) 

a=t,t- m=0,±l 

= [_EEi + ^4E-E]0 +1 + _E^2 + ( e + M a A 2 )0+i + Ap^+1 + V2Ml A ^o 

+ g a ,+l\^a\ 2 + 9n l^m| 2 )0+l +5 , S (l0+l| 2 + I0o| 2 - |0-l| 2 )0+l + ^ 0 ^ 10000 , (7) 

a=t,4- m=0,=tl 


19 













( 8 ) 


ih ~W = + MM 2 T + + £( ^° + ^ m A((f >+1 + 0-i) 

+(£ fl'a,o|'0a| 2 + 9n \ <fim |“) 4>0 + 9s (10+11" + 10-11 2 )0O + 2g s 0j~</>+l</>_l, 

a=t,4- m=0,±l 

+ Mm ^ r ]0-i + + ( e + MaX 2 )^ - \p x 4>-i + V2m A (p 0 

+ ( fi'a,-l|'0a| 2 + 9n \ 4>m V)^-! + 9s (1 4>-l 1“ + 1001“ — 10+11 2 )0-l + S's^+l^O^O) (9) 

a=t,4- m=0,±l 

The damped projected Gross-Pitaevskii equation. The damped projected Gross-Pitaevskii equa¬ 
tion (PGPE) ^ is used to obtain the ground state of atomic-molecular BECs. By neglecting the 
noise term according to the corresponding stochastic PGPE the damped PGPE for the atomic- 
molecular BECs is described as 

d'lpa = V{- l -Ea^ a dt + (/J 0 - H a )ip 0 d.t}, ( 10 ) 

n /c^i 

dtfom = V{- l -Hm$mdt + y~ H m )(j) m dt}, ( 11 ) 

Ft kgT 

where, H a ^ a = (a =|, j,), = ih^f- (m = 0, ±1), T is the final temperature, ks is 

the Boltzmann constant, //„ and fi m are the chemical potential of atom and molecule respectively. 
7 a ( 7 m ) is the growth rate for the ath component (mth component). To obtain the ground state, we 
set the parameter = 0.05 directly. The projection operator V is used to restrict the 

dynamics of atomic-molecular BEC in the coherent region. The initial state of each component is 
generated by sampling the grand canonical ensemble for a free ideal Bose gas with the chemical 
potential /i m 0 = 2/v 0 . The final chemical potential of the noncondensate band are altered to the 
values /i m = 2 /v 
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Figure 1 Experimental setup for creating the spin-orbit coupled atomic-molecular BECs 

of 87 Rb. (a) The experimental geometry. The formation of spin molecular BECs: laser field LI fin¬ 
ishes the spin dependent photoassociation laser field LT performs the A F = —1, Amp = +1 
transiting process of the excited molecules, and laser field L2 induces the excited molecules to 
emit a photon and forms the stable spin-1 molecular BECs. The creation of SOC: the laser fields 
L3 and L4 (L5, L6 and L7) are used to create the SOC of atomic (molecular) BECs along the x 
axis. B denotes a weak magnetic field along the z axis, which is used to induce the Zeeman shift, 
(b) Level diagram of Raman coupling within the atomic BECs. Two lasers (L3 and L4) is used 
to couple states of | t) (i.e., |1, —1)) and | |) (i.e., |1, 0)). This leads to the creation of the SOC 
in atomic BECs. (c) Level diagram of Raman coupling within the F — 1 molecular BECs. The 
three lasers (L5, L6 and L7) couple states of </>_i, 4> 0 and (j) +1 . So the SOC in molecular BECs 
is created, (d) The procedure of creating |F = 1 ,mp = —1) state molecular BEC. A close pair 
of atoms from the condensate in a state | l) absorb a photon from laser field LI, then transit to 
a bound excited molecular state |F = 2 ,m F = —2), i.e., the spin dependent photoassociation 
process |1, —1) <g) 1 1, —1) —> |2, — 2) 21 . Via the laser field LT, the excited molecule transits to 
|F = 1, mp = —1) state. Finally, laser field L2 induces the excited molecular state to emit a 
photon into molecular state 4>-i which keeps F — 1 and rn p = —1. (e) The procedure of creating 
|F = 1, rriF — 0) state molecular BEC. Under laser field LI, the spin dependent photoassociation 
process is 11, —1) <g) 1 1, 0) —> |2, —1). Via laser field LT, the A F = —1, Amp = +1 tran¬ 
siting process leads to excited molecule of |F = 1 ,mp = 0) state. The stable molecular BEC 
(F = 1, mp = 0) is obtained via an emission of photon by laser field L2. (f) The procedure of 
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creating \F — 1, m F = +1) state molecular BEC. The spin dependent photoassociation process is 
11, 0) <£) 11, 0) —> |2, 0). The stable molecular BEC with F = 1 and m F = +1 is obtained via the 
emission of photon by the laser field L2. 

Figure 2 Creating a carbon-dioxide-like Skyrmion in the SOC atomic-molecular BECs. 

(a) The densities of each component when the system reaches the equilibrium state. The com¬ 
ponents are pointed out by the title: | t)> I I), rn F = — 1, m F = 0 and m F = +1, respec¬ 
tively. The particles number (A^, jVj_, N_i, N 0 , N + i) is (2.56 x 10 3 , 2.98 x 10 3 ,1.72 x 10 3 , 3.22 x 
10 3 ,1.59 x 10 3 ). The strength of SOC A = 0.5, the rotation frequency fl ro tation = 0.2a;, where 
ui = 200 x 27 tHz. We set the parameters = g FA = g A with the scattering length a A = 101.8a B , 
. 94 ., 4 . = 0.95 g A , gAM = 0.5 g A , g n = 2 g A , g s = 0.03 g A , x = 0.08, £ = 0, fCt = 0.8 and 
y m = 2 n a = 3.5/i mj0 = 7/i ai0 = 28hu. (b) The corresponding phases of each component in a. (c, 
d) The spin texture of the atomic BECs and molecular BECs, respectively. The color of each ar¬ 
row indicates the magnitude of S z . (e, f) The corresponding spin texture under the transformation: 
(S^S^S’J = (S z ,S y , —S x ). (g) A combined scheme of the carbon-dioxide-like Skyrmion. A 
Skyrmion (center) couples with two half-Skyrmions, which have to arrange at the two sides of the 
Skyrmion. (h) The difference of spin vectors between the atomic BECs and the molecular BECs, 
where y — 0. (i, j) The topological charge density of the atomic BECs and molecular BECs, re¬ 
spectively. The ellipse and rectangles mark out the regions of a Skyrmion and two half-Skyrmions, 
respectively. The symbol o, *, >, * and < are the position of vortices formed by the m F = — 1, 
m F = 0 and m F = +1 components, respectively, (k) The sum of the topological charge densities. 
It denotes the creation of carbon-dioxide-like Skyrmion that the ellipse overlaps with the two rect- 
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angles. The unit of length, strength of x, and SOC are 0.76/zm, hu, u, and 0.96 x 10 3 m/s, 
respectively. 

Figure 3 The effect of SOC on the vortices in rotating atomic-molecular BECs. The 

rotation frequency is Q rot ation = 0.5a;, where u = 200 x 27rHz. We set the parameters = 
9t,l = 9a with the scattering length a A = 101. 8a B , 9u = 0.95 g A , gAM = 0.5 g A , g n = 2gA, 
g s = 0.03 g A , x = 0.08, £ = 0, Cl A = 0.8 and /z m = 2 fi a = 3.5/z mi0 = 7/x a>0 = A2hu. (a) 
Densities of each components under the equilibrium state, the strength of SOC A = 0. Note that 
the components are pointed out by the title: | I m F = —1, rnp = 0 and nip = +1, 
respectively, (b) The corresponding phases of each components in (a), (c, d) The corresponding 
densities and phases with A = 1. (e, f) The corresponding densities and phases with A = 2. (g-i) 
The position of vortices in (a), (c) and (e), respectively. In (g), the rectangle points out a vortex 
cluster. In (h), the rectangle points out the region where the vortices form lattice. In (i), the red 
(blue) rectangle marks out a carbon dioxide vortex, and the black rectangle points out the vortices 
lattice. The particles number (N^, N±, N_i, N 0 , iV +1 ) for the three cases are (5.82 x 10 3 , 8.01 x 
10 3 ,4.96 x 10 3 , 8.23 x 10 3 ,4.93 x 10 3 ), (7.86 x 10 3 , 8.80 x 10 3 , 7.95 x 10 3 , 3.93 x 10 3 , 7.85 x 10 3 ) 
and (11.42 x 10 3 ,12.29 x 10 3 ,10.03 x 10 3 ,2.00 x 10 3 ,10.04 x 10 3 ), respectively. The unit of 
length, strength of x, ^a and SOC are 0.76 /im, hu, u, and 0.96 x 10” 3 m/s, respectively. 

Figure 4 The effect of SOC on the spin texture in the atomic-molecular BECs. (a-c) 
The spin textures of atomic BECs in Figs. 3a, 3c and 3e, respectively, (d-f) The spin textures 
of molecular BECs in Figs. 3a, 3c and 3e, respectively. Note that the spin textures are under the 
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transformation: (S^., S^, S*J = (S~, S y , —S x ). The color of each arrow indicates the magnitude of 
S' z . The blue ellipses point out the regions of the carbon-dioxide-like Skyrmion that Skyrmion of 
atomic BECs couples two half-Skyrmion of molecular BECs. In addition, we mark the position of 
vortices in order to illuminate the relationship between spin texture and position of vortices clearly. 
The meanings of the marks are the same as that in Figs. 3g-3i. (g-i) The corresponding topological 
charge density in Figs. 4a-4c, respectively, (j-l) The corresponding topological charge density in 
Figs. 4d-4f, respectively. In addition, we mark the position of vortices in x < 0 region in order to 
illuminate the relationship between spin texture and position of vortices clearly. The meanings of 
the marks are the same as that in Figs. 3g-3i. (m-o) The total topological charge density of (g, j), 
(h, k) and (i, 1), respectively. We use blue ellipses with the same size to point out the regions of the 
carbon-dioxide-like Skyrmion in (a, d, g, j, m), (b, e, h, k, n) or (c, f, i, 1, o). The unit of length, 
strength of x, ^a and SOC are 0.76 fim, hu, lu, and 0.96 x 10^ 3 m/s, respectively. 

Figure 5 The effect of atom-molecule conversion on the spin textures in the SOC atomic- 
molecular BECs. The rotation frequency is Qrotation = 0.5cu, where to = 200 x 27 tHz. We set 
the parameters A = 1, = 0.8, = g,\ with the scattering length a A = 101. 8a B , 

544 0.95ry\, Qn 25 , a, 5 , s 0.03(jtA, qam 0.5 g /\, & 0 and /r m 2 fi a 3.5 gm,o 

7/Vo — 42 ^cu. (a, d) The spin textures of atomic and molecular BECs with x = 0.0. (b, e) The 
spin textures of atomic and molecular BECs with x — 0.01. (c, f) The spin textures of atomic 
and molecular BECs with x = 0.02. The spin textures are obtained under the transformation 
(S^., S' y , s!.) = (S z , Sy, —S x ). The color of each arrow indicates the magnitude of S' z . In (b, e) 
or (c, f), we use blue ellipses with the same size to point out the regions of the carbon-dioxide- 
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like Skyrmion. [In (a, d), we can not use the ellipses to mark out the region where Skyrmion 
of atomic BECs just couples two half-Skyrmions]. In addition, we mark the position of vortices 
in order to illuminate the relationship between spin texture and position of vortices clearly. The 
meanings of the marks are the same as that in Figs. 3g-3i. The corresponding densities of the 
atomic-molecular BECs are shown in Figs. Sla-Slc [see the Supplementary Information], (g-1) 
The corresponding topological charge densities of the spin texture in Figs. 5a-5f, respectively, (m) 
The total topological charge densities of (g) and (j). (n) The total topological charge density of h 
and k. (o) The total topological charge density of i and 1. In (b, e, h, k, n) or (c, f, i, 1, o), we use 
blue ellipses with the same size to point out the regions of the carbon-dioxide-like Skyrmion. [In 
(a, d, g, j, m), we can not use the ellipses to mark out the region where Skyrmion of atomic BECs 
just couples two half-Skyrmions]. In addition, we only mark the position of vortices in the region 
of x < 0. The unit of length, strength of y, fI a and SOC are 0.76 /im, hu, u, and 0.96 x 10 3 m/s, 
respectively. 
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